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Abstract 

This study is concerned with the decay behaviour of a passive scalar 9 in three-dimensional flows 
having bounded velocity gradients. Given an initially smooth scalar distribution, the decay rate 
d(9 2 )/dt of the scalar variance (9 2 ) is found to be bounded in terms of controlled physical parameters. 
Furthermore, in the zero diffusivity limit, k — > 0, this rate vanishes as K a ° if there exists an ao € (0, 1] 
independent of k such that (\(—A) a / 2 9\ 2 ) < oo for a < ao- This condition is satisfied if in the 
limit k — ► 0, the variance spectrum Q(k) remains steeper than k~ l for large wave numbers k. When 
no such positive ao exists, the scalar field may be said to become virtually singular. A plausible 
scenario consistent with Batchelor's theory is that Q(k) becomes increasingly shallower for smaller k, 
approaching the Batchelor scaling in the limit k — ► 0. For this classical case, the decay rate also 
vanishes, albeit more slowly — like (lnP r ) _1 , where P r is the Prandtl or Schmidt number. Hence, 
diffusion anomaly is ruled out for a broad range of scalar distribution, including power-law spectra 
no shallower than k . The implication is that in order to have a K-independent and non- vanishing 
decay rate, the variance at small scales must necessarily be greater than that allowed by the Batchelor 
spectrum. These results are discussed in the light of existing literature on the asymptotic exponential 
decay (9 2 ) ~ e -7 *, where 7 > is independent of k. 

xxxxxxxxxxxxxxxxxxxxxxxxx 

1 Introduction 

The transport and diffusion of a passive scalar 9(x,t) in incompressible fluid flows u(x,t) is governed 
by the advection-diffusion equation, 

d t 9 + u ■ V<9 = kA9, V • u = 0, (1) 

where k is the diffusivity. The incompressible flows u may be described by a separate system of 
equations, for example the Navier-Stokes system, or just satisfy some prescribed conditions. This 
study considers three-dimensional flows having bounded velocity gradients, i.e. |Vw| < 00. For 
simplicity, periodic boundary conditions are imposed on both u and 9, and the scalar spatial average 
is set to zero, i.e. (9) = 0. The initial scalar distribution is assumed to be smooth. The diffusive decay 
of the scalar variance (9 2 ) is governed by 

±(9 2 ) = -2 K (\V9\ 2 }. (2) 

A problem of long-standing interest is to what extent this decay can be accelerated as the scalar field 
becomes "less smooth" through advective amplification of its gradients. 
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This problem, together with related issues concerning the scalar spectrum, has been actively studied 
for decades, 1-21 most notably since the seminal work of Batchelor and coworkers. 3 ' 4 In one of these 
studies, Batchelor 3 considered presumably non-singular Navier-Stokes flows, where the fluid viscosity v 
is much greater than the scalar diffusivity re. This so-called regime of large Prandtl or Schmidt number 
P r , defined by P r = v/k, is of both practical and theoretical interest, where a fundamental question is 
whether the decay rate d(6 2 )/dt would remain non-zero in the limit re — > 0. The general belief, having 
its root in the phenomenological theory of turbulence, is that this rate becomes independent of re and 
remains non-zero in that limit (even for two-dimensional flows, which are presumably poorer stirrers as 
compared with their three-dimensional counterparts). Thus, the scalar variance is believed to cascade 
to the small scales and undergo anomalous diffusion just as the fluid kinetic energy. However, doubts 
have been cast upon this picture. 11,16 In particular, Sreenivasan 16 suggested that the analogy between 
the energy and scalar variance is rather weak and questioned whether there is a passive scalar cascade 
at all (for a discussion see Chapter 5 of Davidson's book). 22 Recently, Tran and Dritschel 23 showed that 
for general two-dimensional flows having finite enstrophy and power-law spectra of (6 2 ) no shallower 
than A; -1 , the decay rate 2re(|V(9| 2 ) vanishes in the limit re —> 0. Their calculation took full account of 
the suppression of scalar gradients by diffusive effects. 

This study extends the above result to the three-dimensional case, with an improvement in the 
method of analysis by making use of an analytic inequality. For flows having bounded velocity gradi- 
ents, the maximally achievable decay rate 2re(| W| 2 ) is shown to be bounded from above. Furthermore, 
in the limit re — > 0, the derived upper bound vanishes provided that the scalar field remains marginally 
smooth. More precisely, if there exists a positive cxq independent of re such that (|(— A) Q / 2 #| 2 ) < oo 
for a < ao, then re(|V#| 2 } < Cre Q °. Here, ao £ (0,1] and C is a constant, depending on |Vw| and 
(|(— A) Q °/ 2 #| 2 ), both being bounded and independent of re. This constraint covers all scalar distribu- 
tion having variance spectrum Q(k) steeper than k" 1 at the spectral tail, i.e. at large wave numbers 
k. This is an improvement over the result of Tran and Dritschel 23 for the two-dimensional case, where 
power-law spectra were assumed explicitly in their calculations. When no such positive «o exists, 
a plausible scenario consistent with Batchelor 's theory 3 is that Q(k) becomes increasingly shallower 
for smaller re, approaching the Batchelor scaling k~~ x in the limit re — > 0. For this classical case, the 
semi-analytic estimate of Tran and Dritschel can be applied, yielding a vanishing decay rate no slower 
than (lnP r ) _1 . Hence, diffusion anomaly is ruled out for a broad class of scalar distribution, including 
power-law spectra no shallower than k^ 1 . The implication is that a non- vanishing decay rate in the 
limit re — > necessarily requires the variance at small scales to be greater than that allowed by the 
Batchelor spectrum. This is consistent with recent studies 1 ' 7,14,19 suggesting spectra shallower than 
A; -1 . This analysis when applied to the forced case implies that for power-law &(k), the scalar vari- 
ance diverges in the limit re — > 0. Finally, the present results are discussed in the light of existing 
literature 8,13 ' 17 ' 20 on the exponential decay behavior at large times, when the exponential decay rate 
becomes independent of re and remains non-zero. 

2 Basic estimates 

This section briefly recalls a basic calculus inequality and related estimates, both being used in this 
study. In what follows, 6 is bounded and assumed to be sufficiently smooth so that the quantities 
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under consideration are well denned. 

The Fourier representation of 6(x,t) is 

9(x) = J26{k)exp{ik-x}, (3) 
k 

where 9{k) is the Fourier transform of 9(x) and k is the wave vector. In (3), the time variable has been 
suppressed for convenience. For real a, including negative values, the fractional derivative (—A)"/ 2 is 
defined by 

(-A) a / 2 9(x) = k a 9{k) exp{ifc • a;}, (4) 
k 

where k = \k\ is the wave number. For a < 1, the following interpolation-type inequality holds 

(|W| 2 ) < (\A9\ 2 ) il - a)/i2 - a) (\{-A) a/2 9\ 2 ) 1/{2 - a) . (5) 

This inequality reduces to the Cauchy-Schwarz inequality and an identity for a = and a = 1, 
respectively. For a < 1, Eq. (5) can be proved by an elementary method, 24 ' 25 and the proof goes as 
follows. One has 

J k 2 Q(k)dk = J (fe 4 9(A;)) (1 - a)/(2 - a) (fe 2a 0(^)) 1/(2 ~ Q) dk 

, , x (l-a)/(2-a) / , \ l/(2-a) 

<[ k 4 e(k)dk) I / k 2a @(k)dk) , (6) 

where the Holder inequality with the pair of conjugate exponents p = (2 — a)/(l — a) and = 2 — a has 
been used. Eq. (6) is equivalent to (5), thus proving (5). For convenient application in the subsequent 
sections, let us rewrite (5) as {\V9\ 2 )^ 2 -^/^-^ < (\A9\ 2 )(\(-A) a / 2 9\ 2 ) 1 ^ 1 ~ a ^ and re-arrange the 
factors to obtain 

(|A0| 2 ) " (IV^I 2 )"/^-") ' 1 ' 

Equation (7) implies that if there exists an ao > such that (|(— A) a °/ 2 #| 2 } < oo for a < ao, then the 
ratio on its left-hand side vanishes in the limit (|V#| 2 ) — > oo because the right-hand side vanishes in 
this limit. This is the case if @(k) is steeper than k^ 1 for large k. 

For the critical spectrum 0(fc) oc k~ l extending to infinity, no such ao exists. Nevertheless, the 
ratio (|W| 2 ) 2 /(|A#| 2 ) also vanishes for bounded (9 2 ). Indeed, if @(k) = ck~ l for k G [fci,^] (where 
ki k2) and negligible otherwise (steeper than k~ x for k > fo), then c = (9 2 ) / 'ln(A;2 / 'k\) ■ It follows 
that 

<|W| 2 ) 2 _ (cJk'kdkf ^ 



(|A0| 2 ) " cfg&dk " Hk2/hY 



(8) 



The limit hi — > oo necessarily entails ln(/c2/A;i) — > oo if &i remains bounded. Hence (|V#| 2 ) 2 /(|A#| 2 ) — > 
as claimed. Note that the inclusion of the contribution from A; > foj where 0(/c) is steeper than A; -1 , 
does not alter (8) in any significant manner (see Tran and Dritschel 23 for details). 

The above results imply that given a finite field 9, the ratio (| V9\ 2 ) 2 / (\A9\ 2 ) vanishes for a broad 
range of spectral distribution having divergent mean-square gradients, (|V#| 2 ) — > oo. This includes all 
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Q(k) no shallower than for large k. This type of distribution may be called "regular". The type of 
"irregular" distribution that gives rise to a non-zero ratio (| W| 2 ) 2 /(|A#| 2 ) in the limit (|W| 2 ) — > oo 
includes power-law spectra shallower than k^ 1 and non-power-law spectra having a significant fraction 
of (9 2 ) at small scales. This fraction must be greater than that allowed by a k~ l spectrum and non- 
diminishing. A more quantitative sense of this fraction is given in what follows. For a finite (but 
otherwise arbitrarily large) wave number fc*, let = #< + #>, where the large-scale component 9 < and 
small-scale component #> are defined by 

9<= J2 ®( k ) exp{ifc • x} (9) 

k<k, 

and 

6> = 0(k)exp{ik-x}, (10) 

respectively. In terms of these components, the ratio (| V0\ 2 ) 2 / (\A6\ 2 ) is given by 

(|W| 2 ) 2 _ «|W<| 2 } + (|V0>| 2 }) 2 



<|A0| 2 ) (|A0 < | 2 ) + (|A0 > |2) 



f (|ve<l 2 ) , (|ve>l 2 ) 
V<|Ae>| 2 > 1 /2 ^ (|Ae>|2)i/2 



|Ae<| 2 ) 



+ 1 



(11) 



<|Ae>H 

Now, since (|V6><| 2 ) < k 2 {6\) < oo, the limit (|V6*| 2 ) — > oo requires (|V6>>| 2 ) — ► oo, which in turn en- 
tails (|A0>| 2 ) -► oo because (|A6>>| 2 ) > /c 2 (| V6>> | 2 ). Furthermore, (|A(9<| 2 ) < ^(6> 2 ) < oo. It follows 
that both <| V^< | 2 >/<|^> | 2 >V2 o and (| A6»< | 2 )/<| A6>> | 2 > -> 0. Hence, if the ratio (| W| 2 ) 2 /(|A0| 2 } 
remains non-zero as (|V#| 2 ) — > oo, say bounded away from zero by e > 0, then from (11), it is readily 
deduced that 

< (l^l 2 ) 2 (|v^>l 2 ) 2 M9 x 

e " <|A£| 2 ) (|A^>| 2 )- 1 j 

It follows that 

'^WT^ <13> 

where the left-hand side inequality is the familiar Cauchy-schwarz inequality for the small-scale com- 
ponent 0>. Eq. (13) gives an explicit bound for (# 2 } in terms of e. Thus, the variance at small scales 
k > k* is bounded away from zero if the ratio (|V#| 2 ) 2 /(|A6>| 2 ) is to remain non-zero in the limit 
(|V#| 2 ) ^oo. 



3 Main results 

This section derives the main results of this study, including the proof of absence of diffusion anomaly 
for a broad range of scalar distribution mentioned earlier in § 1 . The results are then discussed in 
the light of existing literature, suggesting exponential decay behavior (6 2 ) ~ e~ 7 * in the large-time 
regime, where 7 > is independent of n. Finally, the analysis is extended to include a scalar source. 
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For this case, the scalar variance is found to become unbounded in forced-dissipative equilibrium in 
the asymptotic limit. 

The notion of diffusion anomaly involves subtle issues, which might be worthy of some clarification. 
As the product k(|V#| 2 ) would trivially vanish in the limit k — > if (|V#| 2 ) < oo, both the limits 
k — ► and (|V#| 2 ) — > oo are key features. Furthermore, since (|V#| 2 ) may not diverge in finite times, 
the limit t — > oo may also be involved. As it turns out in the present study, (|V0| 2 ) does not blow up in 
finite times. Hence, the limit t — > oo is indeed another key element of the problem. Now, the question 
is in what manner the limit k — ► should be considered to account for the associated divergence of 
(|W| 2 ) in infinite time. 

This study takes an appropriate approach that answers the above question satisfactorily. Given 
fixed k > 0, let us consider the global maximum of the mean-square scalar gradients (|V6>t| 2 ), where 
0T = 9(x,T) and T denotes the corresponding time upon which this maximum is achieved. In general, 
T depends on n and is expected to diverge in the limit k — > as discussed in the preceding paragraph. 
The present study does not attempt to determine how T depends on k, except that the hypothetical 
scaling T ~ lnP r is examined later on in §3.4, where a related issue is discussed. The behavior of 
k(\V0t\ 2 ) is then monitored against all admissible distribution 9t and its upper bounds are deduced 
in a ^-independent manner. Being derived in this way, these bounds are valid for all k, including the 
limit k — > 0. Evidently, both limits t — > oo and (|W| 2 ) — > oo are being nested within the limit n — > 0. 
Hence, the subsequent constraints on diffusion anomaly deduced from these bounds appropriately take 
into account the divergence of (|V(9| 2 ) in infinite time. In essence, this approach is consistent with 
the double limit lim K _^o h m t— >oo taken in a selective way, in the sense that a dynamical property (here 
being the peak mean-square gradients) is registered during the first limit t — > oo and then monitored 
during the second limit k — ► 0. Evidently, without such a selection, this double limit would yield a 
trivial result because no smooth scalar distribution would survive the limit t — > oo before k tends to 
zero. In passing, it is worth mentioning that the asymptotic regime of small but positive k and large 
but finite time, presumably t S> T, is relevant in the study of exponential decay behavior mentioned 
in the opening paragraph of this section. 

3.1 Production and dissipation of scalar gradients 

The evolution of V6> is governed by 

d t ve + (u ■ v) w = u xve - (ve ■ v)« + kAvo, v ■ u = o, (14) 

where u = V x u is the vorticity. Apart from the "Coriolis" term u x V0, whose sole effect is to rotate 
V# without changing its magnitude, Eq. (14) resembles the three-dimensional vorticity equation in 
several aspects. For example, the gradient V# is advected by the velocity field u just as the vorticity u. 
As another example, the term (V0 • V)u, which is responsible for the production of scalar gradients, is 
an analogue of the vortex stretching term. Despite these apparent similarities, there are fundamental 
differences between "genuine" turbulence and passive scalar advection. First, the passive vector V9 
is irrotational whereas u is solenoidal. Second, 9 is materially conserved while there exists no such 
analogue in the vorticity dynamics. Third, genuine turbulence presumably becomes "more turbulent" 
in the inviscid limit while the limit of zero diffusivity in scalar advection does not make the advecting 
flows more effective stirrers. Among these differences, the third one appears to be most prominent, 
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directly related to the active (and nonlinear) and passive (and linear) dynamical nature of u and W, 
respectively. 

Multiplying (14) by V9 and taking the domain average of the resulting equation yields 

lj t (\ve\ 2 ) = -(v9.(ve.v)u)- K (\Ae\ 2 ) 

<||V^(|W| 2 )- K <|A0| 2 }, (15) 

where the advection and rotation terms identically vanish and the inequality is straightforward. Here 
llVulloo denotes the supremum of Vu and, in this study, is assumed to be bounded. This is the 
maximal degree of smoothness requiring of u. Given an initially smooth 9, the manner in which 
(|V0| 2 ) grows as k is decreased determines the outcome of the decay rate 2k(\V9\ 2 ). This depends 
primarily on how rapidly (|A0| 2 ) grows. Note that even though (|A#| 2 ) cannot be controlled a priori, 
it does good rather than harm, effectively giving rise to the boundedness of (|V#| 2 ) for k > 0. This 
will become apparent in the next subsection. 

A well-known fact which will be exploited later is that for sufficiently small k, (|V6*| 2 ) grows 
exponentially in time for an extended period in the initial stage until diffusive effects become significant. 
For bounded UVu^, it turns out that more rapid growth is not possible, even in the limit k — > 0. This 
is apparent from the inherent linearity of (14) and can be stated more quantitatively in the following 
estimate, being obtained by integrating (15), 

(|V0| 2 ) < (| W | 2 ) exp{2 /"* \\Vu\\ ds} < (| W | 2 } exp{20i}- (16) 

J o 

Here 9q = 6(x,0) and £1 is an upper bound for HVw^. Eq. (16) allows (|V6*| 2 } to grow exponentially 
in time at a rate no greater than 2Q and forbids any possibility of a finite-time divergence. 

3.2 Bounds for the variance decay rate 

For a preliminary estimate of the decay rate, let us consider (15) in the form 

1 d 



_(|A0| 2 ) L (M\ 2 ) 2 K n vol 2A 

^ ( n{ ° 2) - K(|w|2) ) ' (17) 

where (7) with a = has been applied to the second equation. The term Q(9 2 ) is bounded and 
decreases monotonically. The relevant initial data is such that (|V#| 2 ) grows initially. This condition 
will be assumed throughout this section. Recall that t = T denotes the time when (|V#| 2 ) achieves its 
global maximum. This means that at t = T, d(\V9\ 2 ) /dt = 0. Eq. (17) then implies that 

k{\V9 t \ 2 ) < n(9^). (18) 

This provides an upper bound for (half) the peak decay rate. Now, since (|V#| 2 ) < (|V#x| 2 ) for all 
t > and (6%) < (0g), one can deduce that k(\V9\ 2 ) < 0(4) < fi(0g) for all t > 0. 
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Subject to certain conditions on the spectral distribution of 9t, which will be stated in due course, 
refined estimates for k(|V0x| 2 ) can be derived. These conditions are admissible but not necessarily 
realizable, except possibly for the Batchelor spectrum, which has been supported by ample numer- 
ical and theoretical evidence. Hence, the results are conditional and as such should be interpreted 
accordingly. Applying (7) to the second equation of (17) yields 

1 d 2 <|A0| 2 ) / (|(-A)^gp)V(^) 2 \ 

2dt {m } - WW) [ <|Vflp>«/a-) " K(|v 1 } J ( 9) 



(\A8\ 



(o(|(-A) a/2 0| 2 ) 1/(1_a) - sdVei 2 ) 1 ^ 1 -")) . 



(|V0| 2 )V(i-<*) 

In (19), the terms ^(K-A) / 2 ^ 2 ) 1 ^ 1-0 ) and k{\ V0| 2 ) 1/(1 ~ a) effectively represent the relative strength 
of advective production and diffusive suppression of (|V0| 2 ), respectively. Their balance shapes up the 
picture of variance decay. As will be seen presently, this balance is shifted in favor of the latter for 
steeper Q(k). Applying the arguments in the derivation of (18) to (19) yields 

K<|V0 T | 2 ) 1/(1 ~ a) < n(|(-A) a/2 T | 2 ) 1/(1-a) . (20) 

Now, let us suppose that there exists an ao £ (0, 1) independent of k such that 

(|(-A) a / 2 T | 2 ) < oo, for q<q . (21) 

It can be seen that if in the limit k — > 0, 0r(fc) = @(k,T) remains steeper than k^ 1 for large k, such 
an ao would exist. Given (21) and upon replacing a by ao, Eq. (20) becomes 

K(|Wrj 2 ) 1/(1 - ao) < fi(|(-A) ao/2 T | 2 ) 1/(1-ao) , (22) 

or equivalently 

k(|V0 t | 2 ) < n (1 - a ^(\(-A) a ^ 2 9 T \ 2 )K a °. (23) 

This gives an upper bound for k(|V0t| 2 ) in terms of k, among other things, particularly ao- The 
divergence of (|V0t| 2 ) obeys the constraint 

(|V0 T | 2 ) < ^ (1 ^ o) (|(-A) Q() / 2 t | 2 )k( Q() - 1 ). (24) 

A remarkable but not surprising feature of (23) is the explicit dependence on ao, which could be 
termed the degree of smoothness of 0t- As it stands, Eq. (23) reflects the obvious physical connection 
between smoothness and diffusion that smoother scalar fields diffuse more slowly. Note that (23) also 
applies to the border case ao = 1. This case is trivial, simply implying that k{\V6t\ 2 ) scales as k 
if (|V0t| 2 ) < oo. Note also that the conditional bound (23) is an improvement over the result of 
Tran and Dritschel 23 for passive scalar in two dimensions and active potential vorticity in various 
geophysical fluid models. A similar improvement is also possible for their case of enstrophy dissipation 
in two-dimensional Navier-Stokes turbulence, for which the assumption UVu^ < oo (which is not 
justifiable) is replaced at a minimal cost by the natural condition of finite vorticity (which is materially 
conserved). By replacing in (7) and (21) by the vorticity and applying the resulting equations to 
their Eq. (3.2), one can deduce an upper bound for the maximum enstrophy dissipation rate that 
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scales as u a °^ 2 ~ a °\ This scaling converges to zero more slowly than the scaling n a ° . The slower 
convergence is the cost for relaxing the condition HVu^ < oo. 

On physical grounds, Qx(k) can become shallower as k is decreased. This means that condition (21) 
may not be fulfilled, i.e. there may exist no a > independent of k such that (|(— A) a 6x\ 2 ) < oo. If 
this is the case, the scalar field may be said to become virtually singular. This includes the Batchelor 
spectrum @(k) cx fe -1 , which has been believed to be the correct scaling of the so-called viscous- 
convective range, i.e. the range between the energy dissipation scale and the scalar diffusion scale. A 
plausible scenario is that 6t(&) becomes increasingly shallower, approaching the limiting scaling k^ 1 
in the limit k — > 0. For this case, applying (8) to the second equation of (17) yields 




2dV' 1 ' - (|W| 2 ) ^ln(A;i/jfc2) 

<|A6>| 2 > / n(0 2 ) a . 

- <|W| 2 > \Hk K /K) «<l w l>l' 

where k\ and k2 have been replaced by the viscous dissipation wave number k u and the diffusion 
wave number k K , respectively. Following Tran and Dritschel, 23 let us identify the ratio fc 2 /fc 2 with the 
Prandtl or Schmidt number P r . By substituting P r = k\jk 2 into (25), one obtains 

2dV l U ~ (|W| 2 ) \ lnP r 



: <|W| 2 ) < ^pLl _ k( | W | 2 ) . (26) 



It follows that 

i 

K.(\Vf)^\ 2 \ < - 

InP, 



k(|V^ t | 2 ) < (27) 



Eq. (27) implies an extraordinarily weak dependence of k{\V6t\ 2 ) on P r . Nevertheless, k(\V9t\ 2 ) 
vanishes in the limit P r — ► oo. Similar lines of calculation 23 for the case of enstrophy decay in two- 
dimensional Navier-Stokes turbulence have found that the enstrophy dissipation rate scales with the 
Reynolds number Re as (lni?e) -1 / 2 . Again, this relatively slower convergence to zero (as compared 
with (lnP,,)^ 1 ) is the cost for relaxing the condition HVu^ < oo. 

3.3 Discussion 

The results in the preceding subsection cover a broad class of scalar distribution corresponding to no 
anomalous diffusion. This class of regular distribution includes power-law spectra no shallower than 
k~ x . The type of "singular" distribution that would give rise to diffusion anomaly would feature a 
significant fraction of {6 2 ) "cascading" to the vicinity of an ever-increasing diffusion wave number 
k K cx k -1 / 2 as k — ► 0. The decay rate would then be proportional to this fraction. As can be seen from 
the estimates of § 2, this fraction would have to be greater than that allowed by the Batchelor k^ 1 
spectrum. This is consistent with recent results 1,7 ' 14 ' 19 suggesting spectra shallower than k^ 1 . Note 
that for such spectra, virtually all (not just a fraction) of (8 2 ) is at the small scales. 

A number of studies 8 ' 13 ' 17 ' 20 have suggested that in the large-time regime, the variance decays as 
(9 2 ) oc e~ 7 *, where 7 > is independent of k. The exponential decay rate 7 differs from the usual 
decay rate 2k(\V8\ 2 ) and is given by 

1 d 2 _ 2k(\V9\ 2 ) 

7 --W)di {e) -^r- (28) 
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In the presence of diffusion anomaly (the exceptional case of singular distribution mentioned above), 
this large-time regime presumably corresponds to t « T(ko) for some small n = ko, beyond which 
k(\V8t\ 2 ) becomes independent of k. If 7 stays approximately constant for a period of time, then the 
ratio (|V#| 2 )/(# 2 ) is required to remain approximately unchanged, i.e. both (|V#j 2 ) and (9 2 ) decay 
approximately at the same rate 7. In the absence of diffusion anomaly, 7 vanishes in the same manner 
as the usual decay rate for t s=s T. For example, for the Batchelor spectrum, Eqs. (27) and (28) imply 
that the exponential decay rate at t = T, denoted by 77-, satisfies 



Although 7t depends on P r and vanishes in the limit P r — > 00, in principle 7 might become sizable for 
t > J 1 , provided that the ratio (|V#| 2 )/(6> 2 ) has increased by a lnP r -fold since t = T. This requires 
of the variance itself a decrease of at least a lnP r -fold, i.e. (6 2 ) < (#f,)/lnP r , because (|V#| 2 ) also 
decays (at least on average if not monotonically) . Note that an increase of (\V9\ 2 }/(8 2 ) requires that 
(6 2 ) decay relatively more rapidly than (|V#| 2 ). This means that the advective production of scalar 
gradients remains considerably active throughout the regime t > T. 

An interesting problem arises from this analysis is that regardless of the realizability of diffusion 
anomaly, the cumulative decay, up to t = T, defined by 



may remain non-zero and be independent of k. A similar problem on the cumulative decay of enstrophy 
in freely evolving two-dimensional turbulence has been briefly addressed in a recent numerical study, 26 
illustrating the absence of enstrophy dissipation anomaly. There, it was observed that by the time 
T(Re) of peak enstrophy dissipation, the enstrophy loss was significant and, more importantly, roughly 
the same for a range of Reynolds numbers Re. It was also noted that T(Re) increased with Re, 
probably as In Re. In the light of these results, it seems plausible that the cumulative decay, defined 
by (30), remains non-zero in the limit k — > 0. For a quantitative sense of how this might be the 
case, let us consider the following picture, which is in accord with the above observations for the 
decay of two-dimensional enstrophy. In accord with (16), (|V6>| 2 ) is expected to grow exponentially 
for a good part of the time interval [0, T], before diffusive suppression becomes significant toward the 
end. Hence, one may write (|V#| 2 ) ~ (|V#o| 2 ) expjciSlt}, where c\ > is a dimensionless constant. 
Given this growth, the requirement of T for the right-hand side of (30) to be independent of k is 
T ~ c^ 1 ^ -1 In(c2-P r ), where C2 > is another dimensionless constant. Indeed, by substituting these 
into (30), the integral on the right-hand side can be readily evaluated, and the result is oc P r . This 
effectively cancels out the diffusion coefficient in (30). Thus the cumulative decay could be non-zero 
and independent of k in the limit k — > 0. In the absence of diffusion anomaly, one may expect this 
picture to be most plausible for the Batchelor spectrum. The reason is that diffusive suppression of 
scalar gradients becomes relatively stronger (as compared with their advective production) for steeper 
spectra, as suggested by the analysis in § 3.2. This means that for power-law spectra no shallower 
than k~ l , the k~ l spectrum is most favorable for the gradient production and therefore highly likely 
to return the most cumulative decay. 




(30) 
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3.4 The forced case 



This subsection considers the dynamical picture of forced-dissipative equilibrium in the limit k — > 0, 
where the scalar variance is perpetually replenished at the large scales by a scalar source. In the 
presence of such a source, the absence of diffusion anomaly would be catastrophic, in the sense that 
the scalar field would diverge, particularly at the injection scales. It turns out that for power-law 
spectra, {0 2 ) diverges, regardless of whether or not diffusion anomaly is realizable. The classical 
spectrum @(k) oc k~ l is the steepest one that could remain bounded while still satisfying the condition 
{0 2 ) — > oo. Steeper spectra would necessarily grow without bound, i.e. blowup at large scales. 

With the addition of a scalar source / to the right-hand side of the advection-diffusion equation 
in (1), Eqs. (2) and (15) become 

lj t (0 2 ) = -K(\V9\ 2 ) + {9f) (31) 

and 

2dt {m ] - WW) 1 w + — — K{m } ) ' (32) 

respectively. Let us consider a persistent source as described above, i.e. {Of) > 0, and forced- 
dissipative equilibrium dynamics, i.e. k{\X79\ 2 ) = {Of). These can be understood in an average sense 
and are assumed for all P r . Suppose that in the limit k — > 0, (0 2 ) remains bounded. Then for power- 
law scaling, @(k) is steeper than k^ 1 . From the above calculation for the unforced case, the first term 
in the brackets on the right-hand side of (32) vanishes as (|V#| 2 ) — > oo. The second term also vanishes 
for a broad class of scalar sources /, particularly those satisfying (|V/| 2 ) < oo. It follows that for 
{Of) > 0, the assumed forced-dissipative equilibrium k{\V0\ 2 ) = (Of) cannot be achieved. This means 
that such an equilibrium requires {0 2 ) — > oo. The spectrum Q(k) = ck -1 , extending to infinity and 
having non-diminishing c, can be seen as a possible solution, corresponding to a logarithmic divergence 
of the scalar variance toward the small scales. This singular distribution is special and anomalously 
diffused. Other cases of blowup via unbounded spectra steeper than k~ l are also permitted. These may 
or may not be accompanied by diffusion anomaly. Note that non-power-law spectra having significant 
variance at ever-smaller scales and {0 2 ) < oo cannot be ruled out. Similar arguments have been used 
by Tran 27 for the problem of enstrophy dissipation in forced-dissipative two-dimensional turbulence 
(in an essentially unbounded domain). There, the vorticity was found to diverge in the inviscid limit. 



4 Conclusion 

In conclusion, this study has derived upper bounds for the peak decay rate of the variance of a passive 
scalar 6 advected by three-dimensional flows having bounded velocity gradients |Vtt| < oo. These 
bounds are expressible in terms of controlled physical parameters and conditional on the distribution 
of the scalar field 0t = 0(x,T) at the time of peak decay t = T. It has been found that in the limit 
of vanishingly small diffusivity k, the peak decay rate vanishes as K a ° if there exists an ao € (0, 1] 
independent of k such that (|(— A)"/ 2 ^ 2 } < oo for a < Qo- This condition is satisfied if at t = T, the 
variance spectrum 0t(/c) is steeper than k^ 1 at the spectral tail, i.e. at large wave numbers k. Such 
a positive ao is not known to exist as numerical simulations 28 indicate an increasing roughness of the 
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scalar field with increasing P r . If this hypothesis turns out to be false, a plausible scenario consistent 
with Batchelor's theory is that @(k) becomes increasingly shallower for smaller k, approaching the 
Batchelor scaling k^ 1 in the limit k — > 0. For this classical spectrum, the peak decay rate also vanishes, 
albeit more slowly — like (lnP r ) _1 . These results rule out anomalous diffusion for a broad range of 
0x(&), including power-law spectra no shallower than k~ l . The implication is that a non- vanishing 
decay rate in the limit k — ► necessarily requires the variance at small scales to be greater than that 
allowed by the Batchelor spectrum. The absence of anomalous diffusion in the above cases can be 
seen as a consequence of strong diffusive suppression of scalar gradients, which are relatively weakly 
produced by advection. In the presence of a scalar source, (6 2 ) has been found to diverge in the limit 
k — > if @(k) obeys power-law scalings. Hence, spectra steeper than k~ l would necessarily grow 
without bound. The classical A; -1 spectrum, being bounded and accompanied by diffusion anomaly, is 
a possible solution, corresponding to a logarithmic divergence of the variance toward the small scales. 
Interestingly, a diverging scalar field would be necessary for transferring the injected scalar variance 
to the small scales for disposal. 

The present results have been derived on the assumption of bounded velocity gradients. This 
assumption is plausible for the Batchelor turbulence regime, where the Prandtl number P r tends to 
infinity while the Reynolds number Re remains bounded. In fact, there exists neither theoretical 
nor numerical evidence suggesting otherwise. The case where Re — > oo and P r > is interesting and 
presumably much more involved. In this case, one would expect the advecting flows to become increas- 
ingly "rougher", possibly leading to unbounded velocity gradients in the limit Re — > oo. Intuitively, 
diffusion anomaly is more plausible for this case because the scalar gradient production term can no 
longer be controlled a priori. This is no surprise as the original notion of energy dissipation anomaly 
is associated with the limit of infinite Reynolds number. 
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